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Abstract 

Two-dimensional Maxwell- Vlasov equilibria with finite electric fields, ax- 
ial ("toroidal") plasma flow and isotropic pressure are constructed in plane 
geometry by using the quasineutrality condition to express the electrostatic 
potential in terms of the vector potential. Then for Harris-type distribution 
functions, Ampere's equation becomes of Liouville type and can be solved 
analytically. As an example, a periodic "cat-eyes" steady state consisting 
of a row of magnetic islands is presented. The method can be extended to 
(toroidal) axisymmetric equilibria. 
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Equilibrium is the starting point for stability and transport studies of as- 
trophysical and laboratory plasmas. In the framework of collisionless kinetic 
theory, equilibrium states should be constructed as self consistent solutions 
of Vlasov and Maxwell equations. To this end, the knowledge of constants of 
motion for the particles in the continuum approximation (microfluids) is of 
crucial importance because then the general solution of Vlasov equation can 
be written as an arbitrary function of the complete set of constants of motion. 
This is feasible only for one-dimensional equilibria, i.e. in this case the three 
constants of motion are the energy {Hg = {l/2)msv'^ + and the two 

canonical momenta {pxs = 'tIs'^x + Q's^xd/) and = msVz + qsAz{y)); conse- 
quently, the distribution functions are of the form fg = fsiHg.Pxs^Vzs)- Here, 
(x, y, z) are Cartesian coordinates, v"^ = v^+Vy+v"^, $(?/) the electrostatic po- 
tential, Ax{y) and Az{y) the components of the vector potential, and the sub- 
script s denotes the particle species. Unlikely, for two-dimensional equilibria 
the complete set of constants of motion is missing, i.e. only the energy Hg and 
the momentum p^s = msVz + qsAz{x, y) conjugate to the ignorable coordinate 
z are known out of the four constants of motion. A good number of solu- 
tions for one-dimensional [I]- [7] and two-dimensional [5]-[TU] equilibria were 
constructed on the basis of modified Maxwellian distribution functions of the 
forms fs = exp{-PsHs)gs{Pxs,Pzs) and fs = exp{-PsHs)gs{Pzs), respectively, 
with arbitrary functions of the conserved momenta and (3s = ^/{kBTs)- 
These equilibria concern neutral plasmas in connection with a special set 
of distribution functions such that it holds (for one-dimensional equilibria) 
Ni{Ax,Az) = Ne{Ax, Az), where Ni (Ne) is the ion (electron) density [1]; 
viz. in addition to the usual quasineutrality condition it was assumed that 
Ni{Ax,Az) is the same function of A^ and Az as Ne{Ax, Az), thus leading 
to a vanishing electrostatic potential. On physical grounds this additional 
assumption is oversimplifying because it ignores the mass difference of ions 
and electrons. Also, finite electric fields associated with macroscopic plasma 
(ion) flows are important in laboratory fusion plasmas for the transitions 
from low to high confinement modes of operation. 

Aim of the present study is to construct analytically a class of quasineu- 
tral, two dimensional Maxwell- Vlasov equilibria. This is accomplished by 
employing the quasineutrality condition (without the additional assumption 
of functionally identical ion and electron densities) to express $ as a function 
of Az{x, y). A similar method was employed by Mynick and coauthors [11] to 
construct numerically by an iteration algorithm one-dimensional quasineu- 
tral equilibria. Also, the method was reviewed recently in Sec. II of Ref. 
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We consider a plasma of electrons and protons at equilibrium with a 
current density in the axial ( "toroidal" ) z-direction. Consequently, the vector 
potential has a single component Az{x,y). In addition to the "poloidal" 
magnetic field with components Bx{x,y) and By{x,y) associated with A^, 
we include for stabilizing reasons a constant axial magnetic field B^o which 
otherwise does not affect the equilibrium. Furthermore, we employ Harris- 
type distribution functions. 
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where Vzs are constant average (fluid) velocities and tt-qs reference densities 
corresponding to Maxwellian distribution functions {Vzs = 0). Using the 
quasineutrality condition, = Nf., where 
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exp(-gs/5s$) / exp (-/^^m^^Vs) e:!cp{l3sVzsPzs) d^v 

J —oo ^ ' 



(with qi = e and qe = — e), the electrostatic potential can be expressed in 
terms of Az{x, y) as 



^{Az) = log 



exp AzeVze(3e - T^'^eV^Je + AzCVzA + TT^^i^iA 



e(/3e+/3i) 



(2) 



Note that for (1) all the integrations of interest in the velocity space can be 
performed analytically. Using ([1]) and ([2]) one flnds for the current density 
{3z = T.sqsS^ooVzfsd'v): 
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1 , , Pi (AzeVzePe - lm,V,l(3e + AzeVz^(3i + lmiV,^,(3, 
-Vzi {2Aze + miVz^) Pi ^ WTp- 
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Therefore, Ampere's equation (V'^Az{x, y) = —fiajzi^z)) assumes a Liouville- 
type form 

\/'^Az = aexp{bAz) (4) 
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where 



a 



exp 



(5) 



b 



e{y^^-v,e)m 
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The equihbrium for = is static, viz. only the electrons are in non- 
thermal motion to produce the axial current. For 7^ 0, however, there is 
a constant ion-fluid axial velocity, J^^v^fid^v/Ni = Vzi, and the ion motion 
contributes to jz- It may be noted here that for distribution functions of 
the form fs{Hs,Pzs) it is not possible to create poloidal plasma velocities 
because of the two missing constants of motion. Even the third constant of 
motion found in Ref. [13] near the magnetic axis does not help to this end 
because poloidal flows vanish on axis. For the equilibrium constructed here 
the pressure is isotropic, i.e. the pressure tensor. 



is diagonal with P^^x = Pyy = Pzz = P (see also Fig. 3). For Vzi = Vze it fol- 
lows that the current density vanishes {a = b = 0) and (jlj) reduces to Laplace 
equation. Therefore, can not be constant on any closed curve in the [x, y) 
plane without being constant in the region within this curve. Consequently, 
the electrostatic potential $ is constant too in this region because of ([2]) and 
the distribution functions become spatially uniform; hence, one recovers the 
well-known equilibrium solution of the Maxwell- Vlasov equations for which 
all quantities are homogeneous. No "confined solution" is possible in this 
case. Also, it is noted here that for V^i = Vze solution ([H]) below, though 
pertinent to an unbounded plasma, becomes singular {^A^ — > 00). 

Introducing dimensionless quantities = x/ L,rj = y/ L, A^ = Az/{BzoL), 
a = aL/BzQ, b = bBzoL with L a length scahng parameter), the general so- 
lution of (jl]) is given by [TH [T^ 
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where u{^,ri) and v{^,ri) are real conjugate functions resulting from + 
irj) = u + iv, with w a differentiable arbitrary generating function. As an 
example of complete equilibrium construction we consider here the function 



with C = C + ^^7; cind e a parameter such that |e| < 1. For this choice of w, 
(IZj) acquires the form 



The function Az labels the magnetic surfaces. The equilibrium configuration 
shown in Fig. 1 consists of an infinite row of identical periodic islands known 
as "cat-eyes" (see for example Ref. [IS]). The islands have magnetic axes at 
C,a = 2k'7i, rja = and a separatrix with x-points at = (2/c + l)vr, rj^ = 
where k an integer. The ordinates of the separatrix are located at C, = C,a, 
7] = ris = i: arctanh (1 + 2e) (see Fig. 1). The equilibrium has the following 
free parameters: tt-qs, Ps, Vzs (s = e,i), Bzo, e, and L. The dependent 
parameters a and b (Eqs. ([5] and (Q) relate features of the distribution 
function to macroscopic equilibrium characteristics (Eq. ([8])). For Vzi ^ Vze 
the physical quantities (B, E = —{d^/dAz)'VAz, jz and P) are everywhere 
regular and vanish as y tends to infinity except for Ey which in this limit 
approaches a finite value. Profiles of and Ey are shown in Fig. 2 for the 
following fusion relevant values of the free parameters: n^i = rioe = lO^^m"^, 
ksTi = ksTe = IkeV, Vzi = lO^m/sec, Vze = lOVz/, i?^o = ^T, L = Im, and 
e = 0.6. Also, the ^/-profiles of and P have an extremum on the magnetic 
axis (Fig. 3). For e = the configuration becomes one-dimensional; this is as 
an extension of the Harris sheet equilibrium [T7j (usually employed as initial 
state in reconnection studies) with finite E and constant axial velocity. 

Quasineutral equilibria with sheared axial flow which may be more per- 
tinent to the improved conflnement modes can be constructed by the alter- 
native distribution functions 
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with /3s and f3zs constants. A similar procedure then leads to A^-dependent 
average axial velocities: 



Ns 



V 



rrisiPs - Pzs) 



and Ampere's equation assumes the form 



(10) 



where the parameters ai and bi are known functions of nos, 13 s and l3zs- Eq- 
( JTOj) . higher nonlinear than (j4]), should be solved numerically. 

In summary, using the quasineutrality condition to express the electro- 
static potential in terms of the vector potential and Harris-type distribution 
functions (Eq. ([1])) we have constructed a class of plane, two-dimensional 
Maxwell- Vlasov equilibria with finite electric fields, constant axial plasma 
velocity and isotropic pressure. The equilibrium was exemplified by the cat- 
eyes solution. Equilibria with sheared axial flow can be derived by alternative 
distribution functions, e.g. ([9]). The method can also be applied in (toroidal) 
axisymmetric and helically symmetric geometries. 
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Figure captions 



Fig. 1: A^-lines of the cat-eyes solution ([7]) as intersections of the magnetic 
surfaces with the poloidal plane. 

Fig. 2: Profiles of the electric field components Ex and Ey associated with 
the cat-eyes solution ([7]). The profiles E^^y) and Ey{x) have chosen at x/L = 
7t/2 and y/L = 0.5, respectively, because E^i^x = 0,y) = Ey{x,y = 0) = 0. 

Fig. 3: y-profiles at x = of the current density, jz, and the pressure, P, 
associated with the cat-eyes solution ([7]). 
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Figure 1: A^-lines of the cat-eyes solution ([7]) as intersections of the magnetic 
surfaces with the poloidal plane. 
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Figure 2: Profiles of tlie electric field components and Ey associated with 
the cat-eyes solution ([7]). The profiles Ex{y) and Ey{x) have chosen at x/L = 
7i/2 and y/L = 0.5, respectively, because E^^x = 0,y) = Ey{x,y = 0) = 0. 




Figure 3: y-profiles at x = of the current density, j^, and the pressure, P, 
associated with the cat-eyes solution ([7]). 
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